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ON THE NOTION OF EXPLORER XI AROUND ITS CENTER OF MASS 


1 


by 

p 

G . Colombo 


Summar y. — In this paper we evaluate separately the 
magnitudes of the induced and the intrinsic magnetic 
dipoles of Explorer XI (Satellite 1961 Nu) needed to 
explain quantitatively and qualitatively the pro- 
cessional motion of the tumbling axis as an effect 
of both the gravitational torque and the torque that 
the earth ! s magnetic field H exerts on the satellite. 
In addition, a good correlation between the first 
derivative of the frequency and the square of the 
component of H normal to the tumbling axis averaged 
over one orbital period has been found. 


The rigid -body motion of Explorer XI around its center of mass can, 
after the complete damping of the spinning motion, be represented as a 
rotational motion around an axis defined by the vector Q, normal to the 
figure axis Z of the satellite; the Cl - axis is perturbed by several 
effects. The observational material (Naumann, 1961; Naumann, Fields 
and Holland, 1961) clearly indicates that the motion of Q is precessional 

(see fig. 2 ). In addition, within the limit of accuracy of the observa- 
tions, the position of Cl with respect tc the body was observed to be 

^This work was supported by a grant from National Aeronautics and 
Space Administration. 

2 

Consultant, Smithsonian Astrophysical Observatory, Cambridge, Mass, 
on leave of absence from the University of Padua, Italy. 
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Figure 3 -"Period of tumbling motion (sec) versus time 






fixed along the axis of maximum moment of inertia (private communication). 
Finally, a slow decay of the angular momentum L Cl, which is a second-order 
effect, has been precisely detected and studied, showing interesting 
features (see fig. 3 ). We will distinguish the first-order effects 
(motion of Cl) from the second-order effects (slow-down of the tumbling 
mot ion ) . 

The processional motion of Cl can be explained by a torque whose com- 
ponents are quadratic functions of the components of Cl, that is, a torque 
with the same geometrical characteristics as the gravitational torque. 
However, the magnitude of the latter, which can be very well determined, 
is between one -half and one -third of the needed amount; and the orienta- 
tion is quite different* 

In this case of a satellite with a very large perigee distance and 
a complicated shape, but without any large amount of hard magnetic mater- 
ials in its structure, the aerodynamic, eddy-current and hysteresis 
torques are two orders of magnitude less than the gravitational torque 
(Wilson, 1961). Only the magnetic torque coming from the interaction 
between the earth's magnetic field and the equivalent magnetic dipole of 
the satellite appears sufficient to explain the motion of Cl, as we have 
already observed in the case of Explorer IV (Colombo, 1961). 

The torque produced by the earth's magnetic field is 

M = I x I, ( 1 ) 

where I is the vector of the equivalent magnetic dipole of the satellite. 
Neglecting the second-order effects, we can write for the equivalent 
magnetic dipole of the satellite 

I=I 0 +m(H), (2) 
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vhere I is the intrinsic magnetization vector frcm the magnetized can- 
ponents and interior current loops; and m(H) is the magnetic dipole in 
duced hy the inte action of the earth* s magnetic field with the ferro- 
magnetic components. It seems reasonable to suppose that the permeable 
components of the satellite have, in the complex, the same property of 
symmetry with res]3ect to the figure axis Z of the satellite, as the ex- 
terior shape has. Therefore, neglecting hysteresis effects, we may 
write equation (l) as 


M = I x H + 


^(h-k)k 


5 (H-K x Q)K X Q + ^ ( 5 - 3)5 (x H. (1 * ) 


Here jj^, p-g, and are coefficients re Is, ted to the geometrical shape and 
the magnetic property of the permeable ccmponents . 


When we average over a tumbling circle we will find that 


m = (1-3)3 x h + h*(h-H) 3 x 3, 


(i") 


where p,* = [j, - ^(p,^+p,g) . It seems to me important to observe that for 

evaluating p^, p.^, and we must take into account that, while p^ and 
p, are related to a component of the earth*s magnetic field sinusoidally 
changing in a tumbling period (13 seconds), p,^ is related to a very slowly 
changing component of the same field. 


Before the launching no measurements of the intrinsic magnetic 
dipole of the payload were made; only a crude evaluation of the intrinsic 
magnetic dipole of the last stage of the rocket was carried out. There- 
fore, we do not know I with enough accuracy. All we know about the 
permeable structure of the satellite is that the material of the last 
stage of the rocket is 4 10 stainless steel, that small elongated cylinders 
of permalloy material were put into the j)ayload for shielding purposes, 
and finally that there is an iron anular plate in the tail, lying in a 
plane normal to the Z-axis • We do not know anything about the permeable 
structure of the radio-transmitter complex. In order to obtain a possible 
explanation of the observed precessional motion of 0, we must examine the 
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values of I_. and p,* that give the test agreement with the observations . 

For this purpose we need better observations of the orientation of Q and 
a precise numerical integration. Since the probable error in the orienta- 
tion of Q is of the order of a few degrees, we studied separately, using 
the very quick averaging procedure, the cases where 1^ or p,* are negligible. 
In both we found the possibility of having good agreement with the ob- 
served motion of ft, if p* = 4tt x 4.76 x 10 e.m.u., 1^ = 0 and also if 

p* = 0, = 0.63 amp - m 2 . While the needed value of p* seems too high, 

an estimate of the actual value is very difficult. The needed value 
of 1^ in the second case is certainly not in good agreement with the 
measurement made on the last stage of the rocket before launching, although 
no measurements have been made for the payload. Our opinion is that per- 
haps the contribution of both magnetic dipoles , the induced and the in- 
trinsic, are significant. In any case, our goal is to give a very easy 
method for a first-approximation study of the phenomenon. In a second 
stage this method may be improved by numerical integration on a high 
speed computer. 

We assume that the following hypotheses are satisfied: 

(a) the axis of rotation of ft is almost fixed in the body in the meaning 
we will state precisely later; 

(b) the permeable structures of the satellite are such that we can write 
equation (l 1 ) neglecting second-order quantities; 

(c) the aerodynamic, hysteresis, and eddy-current torques are negligible 
for the explanation of the main precessional motion of ft. 

We prefer to leave (a) as an hypothesis since we think a theoretical 
study of the mechanism of the stabilization of the ft-axis in a region 
very close to the axis of maximum moment of inertia would require a 
good knowledge of the internal dissipation of energy (nutation damper) . 

The ft-axis cannot be precisely fixed in the body: if it were, it 

would also be fixed with respect to a fixed reference system and we 
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would, not have the observed precessional motion. However, the total an- 
gular-velocity vector of the satellite may be considered as the sum of 

— ¥ 

the angular velocity vector oo Q of the tumbling motion, fixed in the body, 
and the precessional velocity of the vector 0, with respect to a fixed 
reference system. The magnitude of this second vector is of the order of 
10 y of the value of u). This means that the total velocity vector may 
always be very close to the tumbling axis (less than one second of arc 
separation), which is consistent with the observed motion. 

In hypothesis (b) we prefer to leave the parameters undetermined 
in view of our poor knowledge of the magnetic properties of the permeable 
components and the fact that these components are moving in the weak 
magnetic field of the earth. As we have already stated Dr. C. Lundquist 
made, before launch, a crude evaluation of the magnetization of the 
last stage of the rocket; his results gave a magnetization vector with 
a large component in the direction of the Z-axis and a small component 
(l/20 of the former) in a transversal direction that was not very 
well defined. For hypothesis (c) we can say first that the magnitudes 
of these torques are of two orders of magnitude less than the gravita- 
tional torque, and second that their orientation would not be in good 
agreement for the explanation of the observed variation in the direction 
of Q. The aerodynamic torque can be fairly precisely determined using 
the hypothesis of neutral drag; the hysteresis and eddy-current torques 
cannot he determined so precisely. In any case, however, all these 
torques are dissipative, and the dissipation process involves torques of 
the order of 1 dyne -cm and not of one -hundred dyne -cm; the latter amount 
would be needed to explain the variation of the orientation of Q without 
a dissipation of the same order of magnitude. 



The validity of the following procedures for the deduction of the 
equation of motion is postulated: 

(d) to compute the torque acting on the satellite, we average over one 
period of rotation of the body around fi (tumbling period); 

(e) we average the torque over one orbital period of the satellite; 

(f) finally, we average over one day. 

The procedures followed in (d) and (e) are the usual ones used in 
the perturbation method for determining the gravitational and aerodynamic 
torques (Beletsky, i 960 ). In one tumbling period, the center of mass of 
the satellite will move along a 100-km arc of the orbit . In the case of 
Explorer XI, this amount corresponds to a l°-variation in the orientation 
of the radius vector from the earth's center E to the satellite's center 
of mass G; that means that the variation in the field is two orders of 
magnitude less than the intensity of the field. Furthermore, in one 
orbital period the variation of the orientation of Q is of the order of 
0°5. We prefer to use procedures (d) and (e) as working hypotheses since 
they also give good results for the gravitational torque. Dr. Leland 
Cunningham made for the Huntsville Center a step-by-step integration 
of the original equations to determine the effects of the gravitational 
torque on Explorer XI. His results were the same as those obtained by 
the averaging procedure (private communication from Dr. Lundquist). As 
for procedure (f), we prefer to use the same averaging procedure. There- 
fore, we are able to arrive quickly at the interesting results that follow. 

At the end of this report we will make a preliminary analysis of the 
slow down of the spinning motion. This slow down is definitely a con- 
sequence of the eddy-current torque and hysteresis effects, since the 

aerodynamic torque is very small (less than 0.1 dyne-cm). Both effects 

— * 

are proportional to the square of the component Hj_ of H normal to the 
tumbling axis. The correlation between the square of this component 
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averaged over one orbital period of the satellite and the first deriva- 
tive of the tumbling period states the nature of the breaking torque, 
even if it seems; more complicated, but not hopeless, to distinguish be- 
tween the two effects (Wilson, 1961 ) . 

A more detailed analysis has been made at M.I.T. This analysis of 
the observations, (a dozen per day) of the variation of the period seems 
to show definitely a term with the period of one day, which we think is 
correlated with the variation in one day of the position of the earth’s 
magnetic dipole, which in 2k hours rotates around the earth’s geographi- 
cal axis . 

The averaging procedure we used in our computation is related to the 
accuracy of the observations of 0. In view of the good observations of 
tu, a more accurate knowledge of the value of this parameter may make a 
numerical integration of the general equation of motion worthwhile. 

1 • Magnetic torque from induced magneti c dipole 

Let us compute the effective magnetic torque coming from the induced 
magnetic dipole. From equation (l Tf ), assuming I • ft = 0, we have 


M = p* (H • Q) H x n. 


( 3 ) 


In the usual notation, let 


a(l-e 2 ) 

1-te cos ( (>-ud) 


cos 6 i + sin 9 cos i j + sin 6 sin i 



(*0 


be the vector equation of the motion of the satellite's center of mass 
G. Here r is the unit vector of the direction EG from the earth's center; 
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i is the unit vector of the ascending node in the equatorial plane of the 
earth's equivalent (magnetic) dipole and k is the direction of the axis 
of this dipole. We assume for the earth's magnetic field the usual first 
approximation (Chernosky and Maple, 19 60). 


H = -4 k - 3(k • r)r 


The average value in one orbital period is given by the integral 


m = jp J (h • n)n x h dt 

0 


(5) 


( 6 ) 


Since 


P P — 

d 0 _ 2ra (1-e ) 2 


dt 


Tr 


(7) 


we obtain from equations (3)> (*0* (5)> aad- (6) the following expression: 


M = 7 27 2-. 1/2 

2TTa (1-e ) ' 


2v. 

Hi — - Sx r 2 


J r^( H • 0)H d9. 
0 


(8) 


From equations ( 3 ) and ( k ) we have 


f = n [l-t-e cos ( e-m)] l _ 3 S i n e C os 0 sin i i 


" a 3 (l-e 2 )3 


(9) 


2 • 2 2 

- 3 sin 0 sin i cos i j + (l - 3 sin 0 sin i) k 
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Let Q, , 0 o , te the components of Q with respect to (i,j,k); then from 
-L 2 3 2 

equation (7), neglecting some small terms in e , we finally obtain 


~ ^ 

M " T7 2T972 n x 

a (1-e ) ' 


Q 2 ”♦ ry-f O 2 . ? 

y ciin i q i 4- £1 sin i cos 1 Q j 


§ ±o i 




(l -3 sin^i + ~g sin^ i) k 


( 10 ) 


3 . / , 9.2 


i cos i (l - ^ sin i) (f^k + f^j) > , 


2. Gravitational torque 


We shall nov compute the gravitational torque acting on the satellite. 
Let (If,W,U) he ar. orthogonal reference frame centered at E, where N is 
the unit vector of the orbit's ascending node in the geographic equator, 
and T? is the unit vector of the earth's axis* Also, let i^ be the in- 
clination of the orbit with respect to the geographic equator, and n be 
the unit -vector normal to the orbital plane. Averaging the gravitational 
torque G over one tumbling period and then over one orbital revolution of 
G, we have, finally, 

'V /vi 

5* =■ ^ o) (A-C)(n • n) n X n . (ll) 

Here C is the moment of inertia of the satellite with respect to the 
Q-axis, and A is the moment of inertia with respect to an axis normal to 
Q and 

<«> 

0 r 
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where h is the characteristic constant of the earth's gravitational 
attraction. For Explorer XI, we have 


i Q = 28?8 , 


n = - sin i^ W + cos iQ U ; 


( 13 ) 


and 

2 

p = 3 ~2 (a _ c ) =1.2x10 dyne -cm. ( l 1 *-) 

The components of G with respect to (N,W,U) are the following: 

f G X =1 e n y B (n - nO , 

\ G y = n z + r 6 n x n y ’ <15) 

\ G ,4en o a , 

N z 4 x y 4 x z 

where q , Q , Q are the components of Q with respect to (N,W,U). 
x y z 

3 . Differential equation of perturbed motion 
The equation of motion 

(16) 

dt 
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is now to be projected onto an inertial reference frame and then inte- 
grated with initial conditions corresponding to an observed orientation 
at the chosen initial time. We prefer to project equation (l6) onto the 
moving reference system (N,W,U). We need the equations for passing from 
the reference system (i,j,k) to the reference system (N,W,U). Let I be 
the angle kU; i|r the angle N makes with tae intersection of the geomag- 
netic equator and the geographic equator; and a the angle iN. We have 


first 


ten a = - 


2 tan 1 sin I 


;in ( i n -I ) + tan | sin (I +i ) 


ccs i = cos 1^ cos I + sin i^ sin I cos i|r ; 


and also 


i = cos a N + sin a cos i^ W *f sin a sin i^ U , 

j = (- sin I sin i^ sin a cos t|r - cos I sin i Q sin a)N , 

+ (cos I cos a - sin a sin i^ sin I sin +) W , (l8) 

+ sin I (sin if sin a cos i^ + cos + cos a) U , 

k = sin I sin ♦ N - sin I cos \|r W + cos I U . 


Here I = 11?5; ar.d V = 2rrt rad/day. In our approximation we will obtain 
from equations ( 17 ) and (l8) 


""sin a = -0.41 sin ^ cos a = 1 - 0.08 cos i|r , 

i 

sin i = 0.5 - 0.18 cos cos i = 0.86 + 0.1 cos 


(17’) 
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and 


_ _ _ 

i = (l - 0.08 cos i|r ) N - O .36 sin ifr W - 0.2 sin i|r U , 

/ j = (0.35 sin i(r + 0.02 sin 2\|r)N + (l - 0.08 cos i|r 

( 18 ' ) 

- 0.02 cos 2i|r)W + (-0.04 + 0.2 cos i|r)U, 

\^k = 0.2 sin f N - 0.2 cos ^ W + O .98 U. 

It is necessary now to note that for Explorer XI the reference system 

(N,w,u) is rotating around the U-axis in a uniform motion with an angular 
velocity (regression of the node) of - O.O 87 rad/day. This means that 
equation (ll) projected over the chosen reference system takes the form 



= m (n , n , n , t) + g , 

x x y z x 


= M(n,n,n,t)+G , 

y x y 2 y 


n y t ) + G 
Z Z 


(19) 


The functions G , G , G are quadratic functions of ft , ft , ft with con- 
x 7 y z x y z 

stant coefficients [equation ( 13 ) ] J M^., M , are also quadratic func- 
tions of ft , ft , ft but the coefficients are periodic functions of t 
through i|/, with a period of one day. It would be possible to perform a 
numerical integration corresponding to some initial condition similar to 
the observed conditions of the 35th day after firing. We choose this 
day since after this we have good observations of the precessional motion 
of ft. The amount of work involved in this numerical computation, even if 
worthwhile, suggested to us, that we use a first approximation to compute 
the average values of the components M , M t , M over one day. 
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The observed numerical variation of the orientation of Q is not more 
than 10° per day: in averaging, we consider the orientation of 0 con- 

stant, using the mean orientation for the day. The displacement of Cl 
from the mean value for the day is not greater than 5°* It is difficult 
to evaluate the er:~or involved in this averaging procedure, but in any 
case we think that the approximation is quite good. We prefer to follow 
this method to confirm quickly our feelings about the nature of the 
torque needed to explain the precessional motion. 

4. Evaluation of trie induced dipole 


We shall now discuss and integrate the differential system obtained 
by the procedure explained above . Let us put 


1*4 

Y = ?( ZJW 


( 20 ) 


For our case , a quick evaluation gives 


y = 1.4 x 4rr p* x 10 dyne -cm. 


( 21 ) 


Averaging over one day we obtain 


M 


= v 0 x / 0,'!2 O N + 0.62 0 W + 0 A 7 0 U - 0,25(0 W + 0 U)|, (22) 

1 “ x y z z y J 


It follows that 




2 2 

M = -0.15v Q n + 0«25y (0 - Q ) f 
x y z z y 


m = -o.i5y n 0 + o.25y n n , 

y x z y x 




m = 0.30Y n n - o. 25 y n n 

X y X Z 


(23) 
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To write equation (l 6 ) in the explicit form, we have to evaluate L = A id, 
where cc is the observed angular velocity. We find that 

L = 1.6l6 • 10® gram-cm 2 -sec . (24) 

Choosing the day as the unit of time, now we can write, finally, equation 
(l 6 ) and we obtain 

2 2 

- (o. 15 y + 6 o)fi Q + (o.25y + 52)(q - fi ) -79 Q , 

y z z y y 

- (o.i5y + 90)q n (o.25y + 52)n q + 79 q , ( 25 ) 

x z y x x 

(0.30v + 30)0 Q - (0.25y + 52)0 n ■ 



The differential system in equation (25) has two first integrals; the 
obvious one, 

q 2 + n 2 + fi 2 = i , ( 26 ) 

x y z 

and 

|( 0 . 15 y + 9°)f^ + i(°*30v + 30)o^ + (0.25V + 52)fi^ n z ' 79 & z - E- ( 27 ) 

The intersection of the sphere in equation (26) with the cylinder in 
equation ( 27 ) gives the path of the vertex of the vector fi with respect 
to the rotating reference system (N,W,U). 

In figures 4a and 4b we plotted the projections on the WU-plane and 

— * 

the NTJ-plane of the observed position of the vertex of fi. If we have 


q* = 4 tt x 4.76 10 


5 


m.u. 


(28) 
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Figure k, (a) and (b). — Projection of the observed path of the vertex of 

. 

Q over the xz and yz planes of a rotating reference system (x, 
ascending node, z, earth’s axis). Solid lines are computed 
path, aid dots are observed position. 
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we will obtain 


Y =6.6 X 10 2 ; ( 29 ) 

consequently equation ( 27 ) becomes 

95 C? + 115 Oy - 220 Q y - 79 n z = E. (30) 

Equation ( 30 ) represents a family of hyperbolas. If we choose E = 49-5, 
we will obtain the hyperbola , 'Y , shown in figure 4a. Figure 4b repre- 
sents the projection on the 0 Cl -plane of the path T of 0 corresponding 
to the arc AB of J . The motion of Q with respect to the reference sys- 
tem (N,W,U) is periodic. The good agreement of the observed path emerges 
clearly (see fig. 4c) 

We need also to compare the equation of motion along the path with 

the observations. The easiest way to do this is to compare the observed 

dfi „ . . 

values of — — as functions of fl , with the value of the same derivative 

computed from the third part of equation ( 25 ), which when we take into 

account equation (29) becomes 

dfi . , 

— - = 0.25 Q ft - 0.24 n Q ■ (31) 

dt x y x z 

In figure 4c we plotted the values of the second term of equation ( 31 ) as 

a function of n , evaluated using equations ( 26 ) and ( 27 ) with E =49-5- 
z dQ 

The dots are the observed values of the same derivative — as a function 

Q . The good agreement is evident . 
z 

5 . Evaluation of the intrinsic dipole 

We shall now show that a component of the permanent magnetization 

o 

(evaluated as approximately O. 63 O amp-m ) normal to the Z-axis can explain 
the observed motion of Q. 
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We have again a family of hyperbolas * We choose the constant E in such a 

way that the curve given by equation ( 36 ) passes through the point = 

- 0.8, Q = - 0.6; that is, E = 66.5* 

7 z 

In figure 5a we have the projection of the path of fi on the Q x fi z , 

Q Cl -planes, corresponding to ji* = O.63O amp-m . The same results 

y z 

would be obtained if we assumed the satellite had a residual spin motion 
and a large axial component of the permanent magnetization. 
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Figure 5 


--Same as figure be, corresponding to the hypothesis of a 

const ant component of the permanent magnetization in the 

— * 

direction of Q. 
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While a residual spin velocity of the needed amount was not observed 
(private communication from Prof. Kraushaar of M.I.T.)* the existence of 
the needed component of the residual magnetization has to be postulated. 

6 . Period of tumbling 


In order to explain the variation in the period of the tumbling mo- 

— 4 

tion, we computed the value of the square of the component Hj_ of H normal 
to Q. For a first approximation, we averaged over one day, and then 

H 1 

plotted the value -r- versus time and the value of the first derivative 

If 

of the period in sec/day in figure 6 . The correlation looks good, at 
least for the positions of the maxima and minima. We have to take into 
account both the averaging procedure and the error of the observations of 
Q, which also affects the computed value of 

If we assume the damping torque is the sum of the eddy-current torque 
(which is assumed proportional to the angular velocity cu) and the hystere- 
sis torque (which is considered independent of ou), we can write 


M d =(ocu+u) . (37) 

The magnitude of the torque needed to explain the breaking is of the 
order of 1 dyne -cm. The evaluation of the coefficients a and u is ex- 
tremely difficult and requires good information about the physical 
properties of the conducting and ferromagnetic components of the satellite; 
we therefore plotted only the first derivative of the tumbling period 
versus time. Since the tumbling period varies from 12.4 to 14.6 seconds 
dur ing the one -hundred day interval of observations of the orientation of 

fi. the behavior of 77 cannot be very different, and the position in time 
1 dt 

of the maxima and minima definitely cannot undergo. appreciable changes. 
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A deeper and more detailed analysis of the observed variation of cu, 
even during one day, is strongly suggested by the accuracy of the observa- 
tions and by the interesting correlation of the diurnal periodic term 
with the variation of the position of the earth's magnetic dipole with 
respect to Q. If we take into account that the earth's magnetic dipole 
makes an angle of 11?5 with the earth's geographical axis, the maximum 
displacement of the earth's magnetic dipole will be of 23° in 12 hours, 
which is much larger than the maximum displacement of Q (about 5°) in the 
same period. This fact makes the effect detectable. 

I am indebted to Dr. C. Lundquist and R. Naumann of the Huntsville 
Center and to Prof. W. L. Kraushaar and Prof. G. Clark of M.I.T. for 
the observational material and for fruitful discussion. I am grateful to 
Miss Cara Munf ord for her help in the computation, and to Miss Joan 
Weingarten for her help in the writing of the paper. 
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